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Abstract 

We describe length and time rescaling across the paraelectric to ferroelectric phase transition in 
BaTiOs. Small ferroelectric clusters are unable to grow in the paraelectric phase and exist for a 
short period of time only. The onset of the phase transition invokes a rescaling of both size and 
time lenghts as necessary ingredients for ferroelectricity to come into existence. The growing of 
ferroelectric domains involves correlated chains of Ti-0 dipoles, which causes domain wall rigidity 
and domain shape stability. Ferroelectric nanodomains may host islands of opposite polarization 
(antiferroelectric defects), which appear as an intrinsic feature of the growing BaTiOs ferroelectric 
material. 

PACS numbers: 77.80.-e, 02.70.Ns, 07.05.Tp, 77.84.Dy 
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Ferroelectric materials are broadly employed as capacitors and memory materials. How- 
ever, even for the classical compound BaTiC>3 (BTO), microscopic details on the cubic 
paraelectric (PE) and on the ferroelectric (FE) phases and phase transformations (from cu- 
bic to tetragonal, orthorhombic and rhombohedral phases on lowering temperature) remains 
puzzling. A macroscopic electrical polarization reflects some amount of coherent behavior in 
the atomic displacements. In BaTiOs, the off-center diplacement of Ti atoms creates a net 
electric dipole moment. In terms of domain growth, a critical size |l| may thus be required 
for the ferroelectric state to be stable with respect to the competing electrostatic energy, 
caused by an asymmetric charge distribution. Recently, theoretical [2J and experimental [3|| 
evidence on fundamental size limits in ferroelectricity has beenj^royided, and important in- 

and experiments [zl, [s| 



sights into device scalability have been collected. Simulations 
indicate the existence of dynamic polar clusters already in the paraelectric, cubic phase, as 
local precursors of longer range ferrodomains. This suggests a rescaling of sizes as a fun- 
damental ingredient of ferroelectricity, rather than a size scale within a certain ferroelectric 
state. Rescaling was shown to play a key role in relaxor systems, recently Q|. 

In this Letter, along this line, we perform MD simulations on the phase transition 
(PT) from PE cubic (c) to FE tetragonal (t) in BTO (T c = 393 K) for understanding 
ferroelectric domain formation. For this, we use transition path sampling molecular dy- 
namics (TPSMD)flol] as a means to elucidate mechanistic details of phase transitions 
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12, 



13 . LL4j and to extract scale and time lengths (this work). 



Two main models for ferroelectricity, the displacive UM and the order-disorder model 

QQfl ,,„„„„.., „,.„,,„. 

facts [19j. In the displacive model, the softening of TO phonon modes detected by neutron 
scattering is microscopically connected to Ti leaving the center of the oxygen octahe- 
dra at the transition onset. EXAFS, XANES Q or X-ray diffuse scattering highlight the 
relevance of <111> Ti displacements already in the cubic phase. Molecular dynamics (MD) 
22| indicate a crossover from displacive soft mode to order-disorder dynamics 



simulations 



close to the transition temperature (KNbOs). The implicit time averaging 23J] of NMR ex- 
periments [24I leads to Ti distortion along <100>, whereby the instantaneous displacements 



are along <111> 



21fl. The interaction of off-center Ti displacements and soft modes are the 



25fl, applied to BTO [26fl. Shell model 



main ingredients of the Grishberg and Yacoby model 

1 n 

J| and effective Hamiltonian MD p] point out the prevalent role of chain-like correlations, 
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FIG. 1: Transition path for temperature (a) and field induced (b) PE to FE phase transition. Ti 
displacements Az are shown in red and blue for (up) and (down) domain polarization, respectively, 
(a) 12 3 unit cells box. Clustering is clearly visible in the cubic phase (0-8 ps), just before domain 
growth onset (8-16 ps). (b) Under a static electric field £ the transformation rate is enhanced, 
however certain chains remain antiparallel for a relatively long time (> 25 ps). The box consists of 
9 3 unit cells. 



which are not disappearing above T c . First principles calculations conclude the off-center 



position of Ti in cubic BTO [27||. In a recent DFT-based approach [28|| the fundamental role 
of antiferroelectric (AFE) coupling between Ti-0 chains of opposite polarization is indicated. 

The inclusion of local <111> Ti distortions, coupled to O shifts into short-ranged chains 
is thus a key move for explaining controversial aspects like active Raman signals in the cubic 
phases, x-ray absorption fine structures, phonon modes damping and diffuse scattering. 

While these facts strongly support the central role of domains for ferroelectric properties, 
no clear picture of domain formation across the ferroelectric transition has been provided 
yet. Especially, the role and change of AFE and FE couplings, and the spatial outreach of 
chains in different ferroelectric phases remain unclear. 

To address these questions, we have performed TPSMD simulations on the paraelectric 
to ferroelectric phase transition in BaTiC>3, using the polarizable shell interatomic potential 



29| in the anisotropic NpT ensemble 



of Sepliarsky et al. 

TPSMD is an iterative process and requires an initial trajectory [10]. The latter was 
chosen in a regime of Ti positions shifted off-center from the ideal cubic arrangement, cor- 
responding to the displacive model. This regime is quickly abandoned in the course of 
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FIG. 2: Snapshots of the dynamic structure of BTO in (a) cubic PE phase and (c) tetragonal FE 
phase. The arrows show the relative displacements of Ti atoms from the centers of the oxygen 
octahedra. Note that Ba atoms globally keep their positions and Ti atoms stay on the average in 
the center of a Ba cube. The diffraction patterns on the {111}* reciprocal plane contain diffuse 
scattering in the (b) cubic phase and (d) tetragonal phase (Bragg scattering was subtracted). The 
set of {001}* diffuse planes is clearly extincted in the tetragonal phase. 

simulations, and a clear tendency to grow Ti-0 chains along the c— >t transition appeared. 
In the course of TPSMD a too narrow box would translate into an increased difficulty to 
stabilize the tetragonal configuration on chain growth. We used this intrinsic length resettl- 
ing indicated by the simulation itself to define an appropriate size of the simulation system. 
Successive enlargement of the box dimensions showed that a 12 3 unit cells system is an 
adequate size to study ferrodomain formation. 

The distinction between cubic and tetragonal structures was ensured by an order param- 
eter, based on monitoring cell parameter changes, and on collecting time-averaged statistics 
over Ti displacements (see details below). 

To better cover the diverse (microscopic & macroscopic) aspects of the PE— >FE phase 
transition in BTO, we present in the following a threefold analysis of the simulation results: 
(a) evaluation of instantaneous and time-averaged Ti-displacements (Fig. 2) with respect to 
<111> directions (Fig. 4), for direct comparison with EXAFS measurements (time-resolved 
method), (b) simulation of diffuse scattering (DS, Fig. 2) (c) elucidation of ferrodomain 
growth, domain morphologies and boundaries (Figs. 1 &3), 

Instantaneous and time- averaged Ti displacements. Both PE and FE BTO phases host 
Ti displacements, connected over oxygen into dynamic -(-Ti-0-) m - chains. The Ti envi- 
ronment is essentially rhombohedral in the FE phase, due to the contraction of three Ti-0 
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bonds and elongation of the remaining three (cf. Fig. 2b), a pattern which is locally present 
already in PE BTO. The degree of deviation of Ti displacements from the closest <111> 
direction amounts 12.5 °, in agreement with XAFS results, 11.8 ± 1.1°, and comparable with 



the 10.8° based on the AFE/FE model 



28J . Instantaneous projections of Ti displacements 



indicate <111> shifts in both the PE and FE phases, whereby on time averaging they only 
survive in the t phase, along [001] (Fig. 4). 

Diffuse scattering. The existence of dynamic chains gives rise to the appearance of X- 
ray DS in the form of {100}*-type reciprocal planes intersecting {111}* plane (Fig. 3) [3jJ. 
The extinction of a set of directions (Fig. 2(d)) is due to Ti shift ordering in the FE phase. 
The existence of locally anti-parallel arrangements of Ti shifts highlights the role of AFE 
coupling within both FE and PE BTO. The coexistence of FE and AFE couplings both in 
the FE and the PE phases carries the appearance of DS. A modified pattern due to a changed 
chain timescale and shift ordering in the FE phase is responsible for the extinction of a set 
of {010}* diffuse planes, an effect that can be suppressed upon doping 33]. The nature of 
DS is clearly related in PE and FE BTO. To better understand this point, especially with 
respect to the extent of FE and AFE couplings, we now proceed to a full domain analysis 
across the PT. 

Chains and ferrodomain growth. To trace the progress of phase transition during 
TPSMD, besides cell shape changes, we calculated transverse correlation [sj| and time au- 
tocorrelation of Ti-0 chain formations. We considered time-averaged Ti displacements over 
200 fs. Transverse correlation is an appropriate quantity to monitor the progress of domain 
growth as well as FE and AFE couplings, which appear as vector field features in Fig. l,a. 
Therein, different colors are associated with different Ti displacement directions (Az) and 
overall domain polarization, red for upwards (+) and blue for downwards (— ). The exis- 
tence of polar clusters already in the PE phase (0-8 ps) can be clearly seen, associated with 
AFE/FE interactions. The (fast) time scale of cluster aggregation/decomposition is less 
than 1 ps, for a cluster size up to 2-3 unit cells of locally rhombohedral pattern, supporting 
recent experiments 3, 0]- However, setting in at t >8 ps, a dramatic change in chain dy- 
namics marks the onset of the transition. The dynamically correlated clusters grow larger 
and literally freeze into domains. The overall slow down of the process continues (20-30 
ps) until a smooth landscape of 180° polar domains is reached (t>40 ps). A quantitative 
description of the transition in terms of order parameter changes is given in Fig. 3. Relative 
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variation of lattice parameters amounts to c/a=1.01, in agreement with the experimental 



value of 1.009 



351 ] . The value of spontaneous polarization (Fig. 3) due to the formation 



of domains is (on the average) vanishing, with a hidden tendency of the polarization P z 
component to drift. The non-zero value of the transverse correlation (Fig. 3) along all three 
Cartesian axes already in the PE phase is related to polar cluster formation. The increase of 
transverse correlation, which is a direct indicator of domain growth, surpasses the timescale 
of the abrupt change of dynamics in the system, the latter measured in terms of time au- 
tocorrelation (Fig. 3). The critical slowing down and regime change takes place between 
8-12 ps. At the crossover the system transforms from floppy-dynamic (polar clusters) into 
almost static (domains). During domain formation the z component of the velocity vectors 
v z of Ti positions decouples from the other two and vanishes. In this regime, contrary to 
the PE phase, any further domain evolution in the FE phase requires a whole chain to 
flip. Clearly, there is a marked rescaling of critical lengths in going from PE to FE, already 
noticed above, in terms of size of polar regions, from clusters to domains, as well as with 
respect to range and selection of AFE couplings. A precise size rescaling of chain growth 
correlation lengths and a detailed representation of AFE/FE couplings in the PE and FE 
phases, as well as domain morphologies are the main results of this analysis. Additionally 
our dynamic simulations provide details on the time scale of domain formation and domain 
boundaries (Fig.l, green). In fact, while the PE phase is floppy and may host distant, time- 
decorrelated polar clusters (coarse and mobile domain boundaries), FE domain formation is 
a slower process, with long-ranged correlation lengths within domains and a changed chain 
flip response (narrow and rigid domain walls) . The latter is responsible for trapping narrow 
AFE islands within larger domains (Fig. 1(a), 24, 32, 40 ps). This effect is characteristic of 
the PT, and reflects an intrinsic properties of BTO to accumulate AFE defects within FE 
domains. This is expected to affect the switching response of BTO. 

Chains and ferrodomain growth, £ > 0. To verify this effect we have performed another 
set of TPSMD simulations under an external static electric field £ along z. This was achieved 
by applying different initial values to the two transition branches (very large value for c—*t, 
tiny for t— >c), and incorporating field value reduction (c— >t) or enhancement (t—*c) into 
TPSMD, until an equal value was reached, 60 kV/cm. As the field favorizes t BTO, a 9 3 
unit cells simulation box suffices here. The changes in the order parameters are represented 
in Fig. 3. To a shortened transition time corresponds a smoother change of the lattice 
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FIG. 3: Change of the structural characteristics monitored by different order parameters. Left 
column corresponds to the tajectory of Fig. l,a, right column to Fig. l,b. Colors: blue, green 
and red correspond to the x, y, and z component, respectively Time correlation is the chain auto- 
correlation (0.5 ps). Note that domain growth (transverse correlation [34)]) surpasses the timescale 
of chain mobility freezing. 

parameters. The P z component of the polarization increases up to a value of about 17.5 
/iC/cm 2 , close to the experimental value of 17.0 fiC/cm 2 [35J]. Both transverse correlation 



and time autocorrelation of chain growth indicate a two step mechanism, at 5-15 ps and 
at 15-20 ps. The transverse correlation reveals further domain growth after 20 ps. A pure 
electrostrictive scenario would imply growth of one domain at the expenses of the other 
after 10 ps under the effect of the field. In BTO however, the tendency of placing AF chain 
defects within large FE domains is kept, and survive in the mature, single domain BTO 
after 25 ps, to the extent that the final steps are affected by AF chain dynamics only. 

In conclusion, we have presented an atomistic investigation of the paraelectric to ferro- 
electric phase transition in BTO. Dynamical polar subcritical clusters, already present in the 
cubic paraelectric phase, grow into larger ferroelectric nanodomains while the Ti-0 dipole 
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FIG. 4: Orthogonal projections of Ti shifts (first octant) during phase transition for (a) 0.05 ps and 
(b) 1 ps averaged structures. A clear evidence that different timescales provide different pictures 
of the transition is obtained. 

dynamics is slowing down. Like in relaxors, local polar moments hints at localized phonons. 
Unlike relaxors, size rescaling in BTO restores a lattice dynamics more close to perfect or- 
dered crystals (with the rhombohedral phase as limiting case of perfect dipole order) , while 
in relaxors only the slow down process assists the formation of mesoscopic regions with local 
polarization. Along this line, the size and time rescaling resulting from our TPSMD sim- 
ulations may represent a general scheme for thinking about ordering phenomena in classic 
and relaxor ferroelectrics. 
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